Polymer chains in confined geometries: Massive field theory approach. 
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The massive field theory approach in fixed space dimensions d = 3 is applied to investigate a dilute 
solution of long-flexible polymer chains in a good solvent between two parallel repulsive walls, two 
inert walls and for the mixed case of one inert and one repulsive wall. The well known correspondence 
between the field theoretical 4> 4 0(n)-vector model in the limit n — » and the behavior of long- 
flexible polymer chains in a good solvent is used to calculate the depletion interaction potential and 
the depletion force up to one-loop order. Our investigations include modification of renormalization 
scheme for the case of two inert walls. The obtained results confirm that the depletion interaction 
potential and the resulting depletion force between two repulsive walls are weaker for chains with 
excluded volume interaction (EVI) than for ideal chains, because the EVI effectively reduces the 
depletion effect near the walls. Our results are in qualitative agreement with previous theoretical 
investigations, experimental results and with results of Monte Carlo simulations. 

PACS numbers: 64.60.Fr, 05.70.Jk, 68.35.Rh, 75.40.Cx 

I. INTRODUCTION 

Solutions of long flexible polymer chains in confined geometries such as thin films, porous media or mesoscopic 
particles dissolved in the solution have been extensively studied during last years, including experimental, numerical 
and theoretical investigations. These investigations showed that polymer solutions and binary liquid mixture in 
confined geometries gave rise to a new phenomena not observed in the bulk. The confinement of critical fluctuations 
of the order parameter in a binary liquid mixture leads to an effective long-ranged forces between the confining 
walls or particles immersed in fluid as it was predicted by 1 . It is named critical (or thermodynamic) Casimir force. 
Such fluctuation-induced forces are omnipresent in the nature. For example, such forces arise from the confinement 
of quantum fluctuations of the electromagnetic field and due to the well known quantum-electrodynamic Casimir 
effect 2 . In polymer solutions the reason for this depletion force originates from the presence of depletion zones near 
the confining walls or mesoscopic particles due to the additional amount of entropic energy for polymers confined 
within the slit or between colloidal particles. For entropic reasons the polymer chains avoid the space between 
the walls or two close particles. This leads to an unbalanced pressure from outside which pushes the two walls or 
two colloidal particles towards each other. In the case of addition of the polymer chains to the solvent of colloidal 
solution effective attraction between particles leads to flocculation 3 . Such solvent-mediated flocculation mechanism 
was observed experimentally for silica spheres immersed in the binary liquid mixture of water and 2,6-lutidine . 
Improving of the experimental technique allowed recently even measure with high accuracy the depletion force between 
a wall and a single colloidal particle 7-10 . It should be mentioned, that the case of two parallel walls gives the possibility 
via the Derjaguin approximation 11 to describe the case of big colloidal spherical particle, whose radius R is large than 
the radius of gyration R g and the distance between particle and the wall L. It indicates, that the investigation of the 
case of polymer solutions confined to geometry of two parallel walls is important not only for description of polymer 
solutions confined to film geometry and porous media, but it is also interesting from the point of view of investigation 
of behavior of big colloidal particles in polymer solutions. 

During long period the interaction between polymers and colloidal particles has been modeled by approximating 
the polymer chains as hard spheres 12,13 . But, such approach does not give possibility to describe correctly behavior of 
small colloidal particles in polymer solution and for the case of colloidal particle of the big size the difference between 
theoretical predictions and experimental data are bigger than 10%. In accordance with this more effective were 
approaches which took into account the chain flexibility. For example, for the case of strongly overlapping polymer 
chains as it has place for the case of semidilute solution, the chain flexibility is taken into account via phenomcnological 
scaling theory 14 ' 15 or self-consistent field theory 16 . In the case of dilute polymer solution different polymer chains do 
not overlap and the behavior of such polymer solution can be described by a single polymer chain using the model 
of random walk (for an ideal chain at ^-solvent) or self-avoiding walk (the real polymer chain with excluded volume 
interaction). The last case corresponds to the situation when solvent temperature is above the 0-point (good solvent) 
and polymer coils are less compact than in the case of ideal chains. The remarkable progress in the investigation 
of this task was achieved by 17,18 via using of dimensionally regularized continuum version of the field theory with 
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minimal subtraction of poles in e = 4 — d, where d is dimensionality of the space. But, as it is easy to see 18 , still there 
are a lot of unsolved problems and the question arises : "How to find the theory which allows to explain experimental 
data in a better way?". One of the methods, which up to our knowledge has not yet been applied to this task is the 
massive field theory approach. This method, as it was shown for the case of infinite 19,20 , semi- infinite 21 systems and 
specially for the case of dilute polymer solutions in semi-infinite geometry 22 gives better agreement with experimental 
data and results of the Monte Carlo calculations. In accordance with this, the emphasis of the present work is on 
the investigation of dilute polymer solution confined to geometry of two parallel walls using the massive field theory 
approach in fixed dimension d = 3. 

The most remarkable properties of fluctuation-induced forces is their universality. They are independent of most 
microscopic details and depend only on a few macroscopic properties such as the adsorption properties of the confining 
walls or shape of the particles. In accordance with this we used different combinations of confining walls, i.e. we 
performed calculations for the case of two repulsive walls, two inert walls and mixed case of one repulsive and one 
inert wall. Besides, taking into account the Derjaguin approximation 11 we obtained results for colloidal particles of 
big radius near the wall and compare the obtained results with experimental data 8 . In the case of two repulsive walls 
we found good agreement of our results with results of Monte Carlo simulations 23,24 . 

II. THE MODEL 

We shall assume that the solution of polymer chains is sufficiently dilute, so that interchain interactions and 
overlapping between different chains can be neglected, and it is sufficient to consider the configurations of a single 
chain. Long flexible polymer chains in a good solvent arc perfectly described by a model of self-avoiding walks (SAW) 
on a regular lattice 25 , 26 . Taking into account the polymer-magnet analogy developed by 27 , their scaling properties in 
the limit of an infinite number of steps N may be derived by a formal n — > limit of the field theoretical (p 4 0(n)- 
vector model at its critical point. The average square end -to- end dis tance, the number of configurations with one 
end fixed and with both ends fixed at the distance x = \J (x*a — xb) 2 exhibit the following asymptotic behavior in 
the limit N — > oo 

<R 2 >~N 2 », Z N ~q N N^~\ Z N {x) ~ q N N- {2 ~ a \ (2.1) 

respectively. v 1 7 and a are the universal correlation length, susceptibility and specific heat critical exponents for the 
0{n) vector model in the limit n — > 0, d is the space dimensionality, q is a non universal fugacity. l/N plays a role 
of a critical parameter analogous to the reduced critical temperature in magnetic systems. 

In the case when the polymer solution is in contact with a solid substrate, then the monomers interact with the 
surface. At temperatures, T < T a , the attraction between the monomers and the surface leads to a critical adsorbed 
state, where a finite fraction of the monomers is attached to the wall and form d— 1 dimensional structure. Deviation 
from the adsorption threshold (c oc (T — T a )/T a ) changes sign at the transition between the adsorbed (so- named 
normal transition, c < 0) and the nonadsorbed state (ordinary transition, c > 0) and it plays a role of a second 
critical parameter. The adsorption threshold for long-flexible infinite polymer chains, where l/N — > and c — > is a 
multicritical phenomenon. 

The aim of the present investigations is to describe the behavior of such dilute solution of long-flexible polymer 
chains confined to a slit geometry of two parallel walls located at the distance L one from another in z- direction such 
that the surface of the bottom wall is located at z — and the surface of the upper wall is located at z = L. Each 
of the two surfaces of the system is characterized by a certain surface enhancement constant Ci, where i = 1,2. The 
correspondent effective Landau-Ginzburg Hamiltonian describing such system is: 

L 



+ ^ J d d - 1 rft(v,z = 0) + ^ J d d -\$ 2 {v,z = L), (2.2) 

where </>(x) is an n- vector field with the components <t>i{x), i = 1, n and x = (r, z), ^0 is the "bare mass", vo is 
the bare coupling constant which characterizes the strength of the excluded volume interaction (EVI). The surfaces 
introduce an anisotropy into the problem, and directions parallel and perpendicular to the surfaces are no longer 
equivalent. In accordance with the fact that we have to deal with the slit geometry (x = (r, < z < L)), only parallel 
to surfaces Fourier transforms in d — 1 dimensions take place. The interaction of the polymer chain with the walls is 
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implemented by the different boundary conditions. As it was mentioned above, we consider the case of two repulsive 
walls (the Dirichlet-Dirichlet boundary conditions) 

Ci — > +00, c 2 — > +00 or 0(r, 0) = 0(r, L) = 0, (2.3) 
two inert walls (the Neumann-Neumann boundary conditions) 

d6(r,z). dd>(r,z). „ ,„ 

ci=0, c 2 =0 or VK , z > \ z = = VK Qz J \ Z=L = 0, (2.4) 

and the mixed case of one repulsive and one inert wall ( the Dirichlet-Neumann boundary conditions) 

ci^+oo, c 2 = or 0(r,O) = O, 9 ^*^ \ z =l=0. (2.5) 

The requirement in Eq.(2.4) describing the inert character of the walls corresponds to the fixed point of the so- named 
special transition 21,28 ' 29 in field theoretical treatment. 

In the present case the only relevant lengths are the average end-to-end distance £r = V < R 2 > ~ N v and the 
length L - the distance between two walls. The properties of the system depend on the ratio £/£r. It should be 
mentioned, that the present field-theoretical approach is not able to describe the dimensional crossover from d to 
d — 1-dimensional systems which arises for L « £r. In this case the system is characterized by another critical 
temperature (see, for example, on situation in magnetic or liquid thin films) and moves to a new critical fixed point. 

In accordance with this the present theory is valid for the case L >> £r. Nevertheless, we performed some 
assumptions, which allowed us to describe the region L « £r. 

The well-known arguments of the polymer-magnet analogy 25-27,30 assume the correspondence between the partition 
function Zy(x, x') of polymer chain with ends fixed at x and x' immersed in the volume containing the two parallel 

walls and the two-point correlation function < </>(x)0(x') > in the field theoretical <j) 4 0(n)- vector model at the 
formal limit n — > in the restricted geometry: 



Z||(x,x';iV,L,t;o) =1£^ N (< ^(x^x') > | n=0 ) (2.6) 

Here the r.h.s. denotes the Inverse Laplace transform /j, 2 — > N of the two point correlation function for a system 
modelled via the corresponding Landau- Ginzburg Hamiltonian in the limit, where the number of components n tends 
to zero. N determines the number of monomers of the polymer chain and represents only an auxiliary parameter, the 
trace along the chain and fixes its size globally. The most common parameter in polymer physics to denote the size 
of a polymer chains which can be observable in experiments is i? s 25 , 26 , 30 ): 

R 2 g =X 2 d ^ (2.7) 

where Xd is a universal numerical prefactor which depends on the dimension d of the system. For an ideal polymer 
chains one has \\ = 5 an d f° r three dimensional case N equals -R 2 /2. For the chains with EVI it could be obtained 
within a perturbation expansion 25 . 

The basic element in our calculations is the Gaussian two-point correlation function (or the free propagator) 
< (fi(x)(j)j(x.') > in the mixed pz representation of the form 



G\\{p,z,z';fj, ,c lo ,c 2o ,L) = — ((kq + «o(ci + c 2o ) + ci c 2o )e K ° L - (k 2 - k (ci + c 2o ) + c la c 2a )e K(,L ) 1 

{(k 2 + k (c 1o + c 2o ) + c lo C2 )e K ^ L -\ z - z '^ + (k 2 - k (ci + c 2o ) + c lo c 2o )e- K ^ L -^ z '^ 
+ {n 2 + Ko (c 2o - ci ) - c lo c 2o )e K °( L - z - z ') + (k 2 - k (c 2o - ci ) - c lo c 2o )e- K ^ L - z - z '^), 

(2.8) 

with kq = \J "p 1 -(- /ig, where p is d — 1 dimensional moment. It should be mentioned, that in the case L — > 00 and 
< z, z' « L (or << z, z' < L) the free propagator (2.8) reproduces the free propagator of the semi-infinite model 
(see 21 ). 
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III. THERMODYNAMIC AL DESCRIPTION 



We consider a dilute solution of long-flexible polymer chains with the slit and allow of the polymer coils exchange 
between the slit and a reservoir outside the slit. Thus the polymer solution in the slit is in equilibrium contact with 
an equivalent solution in the reservoir. We follow the thermodynamical description of the problem as given in 18 . The 
free energy of interaction between the walls in such a grand canonical ensemble is defined as the difference of the free 
energy of an ensemble where the wall separation is fixed at finite distance L and that where the walls are separated 
on infinite distance one from another: 

SF = -k B TN In ( Z ' |(L) ) , (3.1) 



Z\\{L -> oo) 

where N is the total amount of polymers in the solution and T is the temperature. Z\\ (L) is the partition function of 
a polymer chain located in volume V containing the walls at a distance L: 

Z\\(L) = j j d d xefVZ||(x,x') , (3.2) 
v v 

with Z\\ (x, x') representing the partition function of a single polymer chain in the slit with its ends fixed at points x 
and x'. For convenience we can renormalise the partition functions Z« (L) and Zn (L — > oo) on the partition function Z 
of one polymer chain in the volume V without the walls. The volume of system V can be divided into two independent 
subsystems Vi and V a which correspond to the volume inside and outside the slit, respectively. This gives possibility 

to split the term ln(— ^ — ) into two parts 

with Z = VZf,, Zt, = J v d x'Zf,(x, x') where Z^ (x, x') is the partition function of one polymer chain in the unbounded 

solution with fixed ends at x and x' and Z Q s{z) = J v d d x' Z\\ (x, x'). 

In the thermodynamical limit (as Af, V — > oo) the contribution from the first term in (3.3) disappear and the 
reduced free energy of interaction Sf per unit area A = 1 of the confining walls may be written as: 

W).,J + J^= Sf ^. 1] , (3.5) 

where n p = Af/V is the number density of the polymer chains in the bulk solution and 

Z HSt (z)= J rf d x'Z ffSl (x,x'), (3.6) 

V„s 

with i — 1,2 and ZnSi ( x ) x ') denoting the corresponding partition functions for a polymer chain in a half space with 
two fixed ends at points x and x'. The functions Zi(z) and ZfjSii z ) depend only on the z-coordinates perpendicular 
to the walls. The reduced free energy of interaction Sf, according to (3.5), is a function of the dimension of a length 
and dividing it by another relevant length scale (namely that for the size of the chain in the bulk, e.g. R x ) yields a 
universal dimensionless scaling function 

&(y)=^f, (3.7) 

where y = L/R x is a dimensionless scaling variable. The resulting depletion force between the two walls induced by 
the polymer solution is denoted as: 
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The total grand canonical free energy £1 of the polymer solution with the slit is: 

tt = -n p k b TALu (3.9) 

with 

L 

u =1 / dz^H. (3.10) 



o 

Taking into account (3.5) and (3.9) we can write for unit surface area A = 1: 

n 



= f b L + f Sl + f S2 +Sf, (3.11) 



n p k B T 

with the reduced bulk free energy per unit volume fb = —l and the reduced surface free energy per unit area 



f Si = J dz (l - . (3.12) 



v HSi 

Further for convenience we can introduce X, the total system susceptibility in the form 



X = 7 / / ddxddx ' < <?i(x)<?i(x') > . (3.13) 



v v 



This definition is consistent with the bulk susceptibility for the unbounded space given as Xf, = -!j to all orders 

of perturbation theory (e.g. 20 ). Zf, being the Inverse Laplace transform of Xb and Zb = 1 to all orders as well. 
Accordingly to (2.6) and (3.13) we can rewrite (3.5) in the form 



Sf =l£^ R2 j 2 {L(X b -X ll ) - Ti - T 2 }, (3.14) 

where X\\ denotes the total susceptibility for a slit geometry and Y$ with i = 1, 2 give two half-space (HS) contributions 
such that / Si = 2X^2^^2/2! Yi} (see Appendix A). 



IV. CORRELATION FUNCTIONS AND RENORMALIZATION CONDITIONS 



Correlation functions which involve N' fields 0(xj) at distinct points Xj(l < z < N') in the bulk, Mi fields 
(j>i(rj 1 ,z = 0) = i^gj^^jj) at distinct points on the wall z = and M 2 fields 2 (rj 2 ,2: = L) = 4> S2 {^j 2 ) & t distinct 
points on the wall z = L, and I insertion of the bulk operator |0 2 (Xfe) at points with 1 < k < I, I\ insertions of 
the surface operator 5^ (R^) at points R^ with l<h<I\ and J 2 insertions of the surface operator i0 2 2 (Rj 2 ) at 
points R; 2 with 1 < l 2 < I2, have the form 21,29 



G (iV',M 1 ,M 2 ,7,7 1> / 2)({x . }){r . i}){r . 2}){Xfe}){R;i}){Ri2}) = 
N' Mi M 2 / h I 2 

< n^( x o n n n 2^ 2(Xfc) n n ^i^) > • w 

i=l ji=l J2=l fc=l li=l i 2 =l 

Here, the symbol < ... > denotes averaging with Hamiltonian (2.2). The free propagator of a Gaussian chain in slit 
geometry in the mixed p, z representation has the form (2.8), as was mentioned above. 

Taking into account that surface fields <f> ai ( T ji) an d surface operators \4> 2 Si (Rj) with i = 1,2 scale with scaling 
dimensions that are different from those of their bulk analogs <f>(x-j) and i0 2 (Xj) (see 21 ), the renormalized correlation 
functions involving N' bulk fields and Mi surface fields on the wall z = and M 2 surface fields on the wall z = L, I 
bulk operators, I\ and 7 2 surface operators can be written as 

r ,(N',M 1 M 2 ,I,I 1 ,I 2 ), rA 
7 -(N' + M 1 +M 2 )/2 7 -M 1 /2 7 -M 2 /2 7 I 7 I ± ? I 2 r ,(JV',M 1 ,M 2 ,/,/ 1 ,7 2 ) / n / 4 o\ 
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where Z^, Z\, Z2 and Z^2, Z^ , Z^ are correspondent UV- finite (for d < 4) renormalization factors. The typical 

bulk and surface short-distance singularities of the correlation functions G^ N > M i> M 2) can be removed via mass shift 
/iQ = /i 2 + Sfi 2 and surface-enhancement shifts Ci a = c% + Sci, respectively 21 . The renormalizations of the mass /x, 
coupling constant v and the renormalization factor Z^ are defined by standard normalization conditions of the infinite- 
volume theory 20 ' 31-34 . In order to adsorb uv singularities located in the vicinity of the surfaces, a surface-enhancement 
shifts Sa are required. In this connection the new normalization conditions should be introduced. It is obvious, that 
in the limit L — > 00 we should have 



lim 



lim 



Gfl' 2,0) (p;fi,v,c 1 ,c 2 ,L)\ p=0 



1 



/=<( .0,2)/ r x. 



p=0 



M + Ci 

1 



M + c 2 

For the renormalization factors Zj,2 where i = 1, 2 we obtain, respectively 



(4.3) 



lim 

L— >oo 



lim 

L— >oo 



9 a(0,2,0)/ r x| 
9 -=v(0,0,2), n | 

~q^2 R Hp;M,«,ci,c 2 ,L)| p=0 



and 



2 ii (/x + ci) 2 ' 
1 

>(Ai + c 2 ) 2 ' 



(4.4) 



lim 

-L — >oo 

lim 

L— >oo 



A(0,2,0;0,l,0) 



G 



(0,0,2;0,0,i; 



(p,P;/i,u,Ci,C2,L) 
(p,P;/i,u,Ci,c 2 ,L) 



p,P=0 



p,P=0 



1 



(^ + ci) 2 ' 
1 

(^ + c 2 ) 2 ' 



(4.5) 



In the limit L — ► 00 all these conditions yield exactly the same shifts <fcj and renormalization factors as in the semi- 
infinite case. It is intuitively clear that in the case of two inert walls or mixed walls situated on big, but finite distance 
L with L > R g such that the chain is still not deformed too much from its original size in the bulk, the shift of 
Cq P — ► c sp may depend on the presence of the other surface and hence on the size of the slit. So, in the case of L > R g 
(or \xL > 1) from (2.8) and (4.3) we obtain new conditions 

1 



lim 
lim 

Lfi>l 



G^' 2 ' 0) (p; /x, v, ci, c 2 , £)| p =o 



G^'°' 2) (p;/x,i;,ci,c 2 ,L)|p = o 



fl + Ci 

1 

M + c 2 



1 + 



2/i zx - c 2 

( 

/i + ci /i + c 2 

2/i /i — Ci 

( 

/i + c 2 /i + ci 



-2/ii 



-2/iL 



+ o( c - 4 ^ L ) 



+ 0(e^ L ) 



(4.6) 



The above mentioned conditions (4.6) give one-loop order corrections to the respective surface-enhancement shifts 
Sci of semi-infinite theory in the case of large, but finite wall separation L . In accordance with this for the case of 
mixed walls we obtain 



<5cf-° = 5 Cl + A< s -°) 



with corrections of order 0(e 2tlL ) 



A (s -°^ = ^(ji + Ce + In 8 - 3 + ln/iL - c 4 " L Ei(-4/iL)^ c -2 ^. 

In the case when both walls are inert, the modified surface enhancement shifts are 

Scf- S = S Ci + A< s " s ) 

with 

-i\ih 



A (s-s) = _ A (s-o) _ ^ ln2 _ 



(4.7) 
(4.8) 
(4.9) 
(4.10) 



The above mentioned corrections 5ci are UV singular for d = 3 dimensions. They provide the singular parts of the 
counterterms that cancel the UV singularities of correspondent correlation functions by analogy as it took place for 
semi-infinite systems (see 21 ). The above mentioned corrections A( s ~°' and A( S-S ) are finite in d < 4 dimensions. 



7 



V. RESULTS FOR GAUSSIAN CHAINS 



Let us consider at the beginning the Gaussian model for ideal polymer chains (vo = 0). As mentioned above it 
corresponds to the situation of a polymer chain under O-solvent conditions. 

For general case of arbitrary ci and c 2 on the confining walls we obtain for the reduced free energy of interaction: 



Sf 



TC 



" A Jr 7 ' 



(p + ci)(/i + c 2 )c^ L - ( M - Cl )( M -c 2 )c-^ 



{4cic 2 - (m(ci +c 2 ) + 2cic 2 )c^ l + (//(ci + c 2 ) - 2cic 2 )c " L } + 



ci 



C2 



/U° \ fj, + c\ fi + c 2 



(5.1) 



First, consider the case of the Dirichlct-Dirichlct (D-D) boundary conditions (2.3) on the confining surfaces. Taking 
the limits — — + oo, — — > oo yields: 



Q D ' D (y) 



(5.2) 



where t = /i 2 L 2 and y — The result indicates that if both being positive, the depletion interaction potential is 
negative and hence the walls attract each other due to the depletion zones near repulsive walls. The inverse Laplace 
transform can only be performed numerically (the plot is shown in Figure 1) or may be expanded for asymptotic values 
of \fr. The obtained results for ideal polymer chains in slit of two repulsive walls are in agreement with previous 
theoretical results obtained in Ref. 18 . But, it should be mentioned, that on plotting these functions the authors of 18 
used a rescaled variable V2R X , which was not mentioned there. 

Now we proceed to the case of two inert walls, what corresponds to the Neumann-Neumann (N-N) boundary 
conditions (2.4). For the free energy of interaction we obtain 



N > N (y) 



0. 



(5.3) 



This corresponds to the fact that ideal chains do not loose free energy inside the slit in comparison to the free chains 
in unrestricted space. The entropy loss is fully regained by the surface interactions provided by the two walls. 

Taking the limits ^ — » oo, ^ — > in accordance with (2.5) (the Dirichlet-Neumann (D-N) boundary conditions) 
from (5.1) we obtain: 



Q u > N (y) = -2y l£ T ^ {2y 2 



1 



1 



r 3/2 1 + e 27? 



(5.4) 



This result can only be evaluated numerically and is plotted in Figure 2. Lets consider different asymptotic regions 
of y. 

Wide slits (y » 1): 

In the case \iL » 1 from (5.2) we obtain for two repulsive walls: 



Q D - D {y)*4y 

The force (3.8) becomes 



erfc f ^ 



1 



y V 7T 



exp 



' 2 



erfc 



(V2y) 



1 



y V27T 



exp (-2y 2 



T u ' u (y)m -4 erfc 



V2 



8 erfc 



And for one repulsive and one inert wall we have: 

e D ^(y)«4 2/ crfc(V2y) 

which implies 



exp (-2y ) , 



r^ Ar (y)«-4crfc(V2y). 



(5.5) 



(5.6) 



(5.7) 



(5.8) 
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These approximating functions are presented on Figures 1 and 2, respectively. 
Narrow slits (y -C 1) 

In the case of narrow slit \iL « 1 the asymptotic solution for (5.2) reads: 



Q D ' D {y)*- 4= + V- (5-9) 



2tt 



and the force simply becomes r D,D (y)xi — 1. 

For the depletion interaction potential (5.4) we get: 



e^G,)* - -4= + y ■ (5-io) 



. 2tt 

For the force we have again T D ' N (y)~ — 1. 

These results can be understood phenomenologically. In our units the quantities and T are normalized to the 
overall polymer density n p . So, the above results simply indicate that the force is entirely induced by the free chains 
surrounding the slit or in other words by the full bulk osmotic pressure from the outside of the slit. No chain has 
remained in the slit. It is reasonable in the case of repulsive walls in the limit of narrow slits. Unfortunately, the 
narrow slit regime is beyond the validity of our approach in the presence of EVI, as mentioned above. But, the 
above mentioned arguments can be used in order to obtain the leading contributions to the depletion effect as y — ► 0. 
We can state that in the case of very narrow slits the chains would pay a very high entropy to stay in the slit or 
even enter it. It is due to the fact that the phase space containing all possible conformations is essentially reduced 
by the squeezing confinement to a size times its original size (for an unconfmed chain). Therefore, the ratio of 
partition function of polymer chain in slit and free chain partition function vanishes strongly as y — > 0, which implies 
directly the function u) in (3.10). Setting uo = and using only the corresponding surface contributions and the bulk 
contribution (/& = —1) in (3.11) must lead to the same asymptotic limits in the narrow slit regime. The advantage of 
this procedure is that no expansion necessary and it should be equally valid in the EVI-regime. 

In Figures 1,2 and 3 the depletion interaction potential Q(y) and depletion force T(y) are plotted for all boundary 
conditions. As expected, the results for mixed walls are located in between the results of two inert walls and those of 
two repulsive walls. 

VI. RESULTS FOR GOOD SOLVENT 

In good solvent the EVI between chain monomers play a crucial role so that the polymer coils occupy the bigger 
volume and are less compact than in the case of ideal polymer chains. The influence of EVI on the depletion functions 
can be obtained in the framework of the massive field theory approach in fixed dimensions d = 3 up to one-loop order 
expansion of the two-point correlation functions G^ 2 ' 0,0 ' restricted in slit geometry (2.2). The bare total susceptibility 
Aj are (see(3.14)) for the slit geometry in accordance with (3. 13), (3. 10) and (4.1) is : 

L L 

X\\ are (Ho,vo,c lo ,c 2o ,L) = i J J dzdz'j(5||(p = 0, z, z'; c 40 , L) 

o o 

L 

v ° I dz " I < ^" ( ~ P = ' z ' z '' ; ^ ' C4o,jL ^ii^ q ' z ''' z '' ;Mo ' Cto ' L ^ii^ p = °> 2; ",z';/4),Cio> i )} • 

q 



n + 2 
6~ 



The two HS contributions denoted by Tj (see (3.14))can be obtained in accordance with (3.12) similarly to (6.1) with 
the propagators of semi-infinite system. Some details for the calculation of these quantities for zero-loop and one-loop 
order for different surface critical points of interest (ordinary, special) are presented in the Appendix A. 



A. Two repulsive walls 



Lets consider first the case of D-D boundary conditions (2.3) on each of the two surfaces. In this case no surface 
divergences appear in the calculation of the correlation functions and any surface renormalization is not necessary at 
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all. Each surface term (f Si , i = 1,2) contributes: 




After performing the standard mass and coupling constant renormalization and additive subtraction at zero 
momentum all divergent terms disappear and the corespondent function XP^ n can be obtained. In order to be 

concise, we do not present here the complicated form for X^en, an< ^ J us ^ discuss the limiting cases of wide and narrow 

slit regimes. 

Wide slits (y > 1): 




FIG. 1: The functions O(y) and T(y) for two repulsive walls 



The massive field-theory approach at fixed dimensions d = 3 gives a rather simple result in one-loop order than 
results obtained in 18 with help of dimcnsionally regularized field theory with minimal subtraction of poles in e- 
expansion. It should be mentioned, that in 18 a wide slit approximation was carried out as well up to the first non 
trivial order (apparently O (e~^ L )). Therefore, we performed calculations up to the next order term <~ O (c~ 2fJ,L ). 
The renormalized total susceptibility for the slit geometry up to one-loop order in d = 3 for polymer case n —> in 
the wide slits regime fiL >> 1 is: 



v D < D r~ 



-1y,L ( g _ Q _ 2 l n : 



In(/iL) 



l n : 



-fj,L 



+ 



4-ln- 



+ e^ L Ei(-/xL) - c^ L Ei(-3/xL) 



•Ei(-4/iL) \. (6.3) 



The exponential integral functions, denoted by Ei(x), can be expanded for large, negative arguments as well in 
accordance with (see e.g. 35 ): e x Ei(— x) = — - + O Thus, for the depletion interaction potential we obtain: 



e(»)w 4 - In 



163 , 3 C E 

ln 

24 2 2 



yerfc 
2 



2y erfc(\/2y) - ^ 



y 

V2 



exp (-2y 2 ) 



2 fy 2 

exp 



7 XC r^ 

4 2^ 



55 
4% 

lnr 

7m 



fc (V2y 



(6.4) 



A comparison of the obtained results to the ideal chain results in a wide slit regime (see Figure 1) shows that the EVI 
reduces the depletion effects for two repulsive walls. 

Narrow slits (y <C 1): 

Following the simple argument obtained from the discussion of the exactly solvable ideal chain model the entire slit 
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contribution u (3.10) to the reduced free energy of interaction 6f in (3.11) is simply set to zero and the depletion 
effect is only calculated from the bulk and surface contributions. In this limit the depletion potential becomes: 



w«v--*r(i- Y 



(6.5) 



and the force again is unity. In Figure 1 one can follow how the two regimes come to match in the crossover regime 
The lowest order expansion in case of wide slits would not be able to show this matching. With these two 
approximations we are in the position to present a rather complete picture of the problem in comparison to the 
approach given in 18 . 

B. One repulsive / One inert wall 

This case has not been studied so far in any approach. Since we are now dealing with an inert wall, the surface 
renormalization should be taken into account. Again, the full result for the renormalized total susceptibility in a slit 
system X\ \ ren has complicated form and we discuss just limiting cases of wide and narrow slits. 

The surface contribution for a repulsive wall coincides with (6.2) and for inert wall we have: 



(6.6) 



Wide slits (y > 1): 

For the total susceptibility up to O (c~ 2aiL ) order we obtain: 



os L 1 ( 2 ln ! ~ 1 

X \\ren L ~-^ ~ ^ ( 1+ 4 



+ 



^ L /21n4 - ln3 - 1 



c- 2 ^ f 31 - 2C E In 3 7 1n2 3 ln(/zL) 
+ ^ 3 I 8 2 4 4/xZ 



. + — Ei(- M L) + — Ei(-4/xL) - — Ei(-3 M L) ^6.7) 



In comparison to the result for ideal chains (5.4) where the lowest order term, contributing to the total susceptibility 
in the wide slit limit is of order O (e _2Ali ), now the additional term of order C(e _Aii ) appears. 



0(y) for 1 inert / 1 repulsive wall 



r(y) with EVI for 1 inert / 1 repulsive wall 
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FIG. 2: The functions &(y) and T(y) for one repulsive and one inert wall with and without excluded volume interactions (EVI) 
In this case the depletion interaction potential becomes : 

(.»! - l) (»e rf c (JL) - oxp (-£)) + ^«fc(VS») 

+ i - c,i - lnl152 ) ( 29erfc i^") - vf oxp <" 2! ' 2) ) + ! I£ -* 



lnr 

^372 



(6.8) 
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The Figure 2 presents the depletion interaction potential O(y) and the force T(y). It clearly indicates that in 
comparison to ideal chains the depletion effect is stronger in the regime of wide slits. 

Narrow slits (y' -C 1): 

Following again the thermodynamic argument, uo is set to zero and only bulk and surface contributions are taken into 
account in (3.11). One gets: 

e (y) ^- (i + ^i- 1 )/!' ^ 

which is also slightly below the depletion potential in comparison to the case of ideal chains (see Figure 2). The 
depletion force is unity. 

Both approximations for wide, as well as for narrow slits suggest the depletion effect to be stronger in the case of 
excluded volume interactions than for ideal polymer chains (see Figure 2). 



C. Two inert walls 



In order to obtain the renormalizcd total susceptibility for a system confined by two parallel inert walls we have to 
apply the surface renormalization scheme suggested by 21 for both surfaces at their surface critical point c 8 ?. Starting 
from (6.1) we obtain for the renormalized total susceptibility: 

f^^- ^3 ( ln2 - \ - ln(l-e- 2 ^)) • (6.10) 

The surface contribution has already been presented in (6.6). Lets consider the asymptotic expansion for wide slits 
[iL >> 1. Taking into account the surface (6.6) and the bulk contributions, the result for the depletion interaction 
potential becomes: 




(6.11) 



This function and its derivative for the force are plotted in Figure 3. 
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FIG. 3: The functions Q(y) and T(y) for two inert walls with and without excluded volume interactions (EVI). Here we 
introduced notations: 8c S emi = Set and 5cf~ s = 8c s iit with i = 1, 2 



It is obvious that only the wide slit approximation can be applied here since the usual argument for the narrow slit 
approximation is no more valid and u> does not necessarily vanish. 

Interestingly, the depletion force turns out to be positive and the walls are repelled from each other. This means that 
the chains rather like to stay in between the slit than leave it. This in turn means that the chains gain enough energy 
from attractive interactions on the walls, which forces them to exert their loss of entropy (due to the confinement) 
onto the walls instead of leaving the slit. 
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It is very instructive to have a more general look on the terms appearing in the free energy of interaction. If now we 
take into account the new normalization conditions for surface-enhancement constants for slit geometry (see Eqs.(4.6)- 
(4.10)), which assume that we have big, but finite wall separation L, the Sf can be written as: 

Here 8c s ~ s is the surface-enhancement constant shifts for slit geometry which appears in the case of finite walls 
separation and Sci is surface-enhancement constant shift in the case of infinite walls separation. In the presented 
approach the same renormalization of critical values Cq P was used and equally the same shifts to the renormalized 
values were obtained. So the first term on the r.h.s. just disappeared on the assumption that the surface-enhancement 
constant shift on one surface is not affected by the presence of the second one. 

In fact this assumption could be doubted and an additional shift through the influence of the second wall (coupling 
effect between the two walls) may appear. Since the interaction potential itself is purely local, such a coupling effect 
can only be mediated through the chain conformations. As a result the number density of monomers near the walls 
might differ in comparison to a semi-infinite constraint and also the shift of the critical point (due to excluded volume 
interactions) can change. As already proposed in 36 this in turn would require a different renormalization scheme for 
the surface critical point, where this coupling effect is to be taken into account. The results of calculations for a slightly 
modified surface renormalization scheme which takes into account the finite surface separation L are introduced at 
the Appendix B and are presented in Figure 3 as well. 

VII. COMPARISON TO PREVIOUS WORK 
A. Theoretical approach 

As was mentioned in the Introduction, the remarkable progress in the investigation of the influence of EVI on the 
depletion interaction and depletion force between two repulsive walls was achieved by ,18 via using of dimcnsionally 
regularized continuum version of the field theory with minimal subtraction of poles in e = 4 — d, where d is dimen- 
sionality of the space. Figure 4 presents comparison of our results obtained in the framework of massive field theory 
at fixed dimensions d = 3 for the case of two repulsive walls and results obtained in 18 . 



Compare ©(y) for 2 repulsive walls with SHKD Compare r(y) for 2 repulsive walls with SHKD 

I 1 . 1 1 1 I 1 1 . . . 




y = LVR X y = L7R X 



FIG. 4: The functions Q(y) and T(y) for two repulsive walls in comparison to 



The results obtained in the framework of both analytical methods are in quantitative agreement. But, one notes 
that the reduction of the depletion effect due to excluded volume interactions is less stronger within the massive field 
approach as compared to an e-cxpansion in one-loop order. It should be noted, that we extended our results up to 
the next e~ 2Ali order. This allowed us to obtain good matching with approximating results in narrow slit limit (see 
Figures 1,2 for 6{y)). 
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B. Simulations 



One of the possibilities to test reliability of the obtained analytical results is to compare them to results obtained by 
Monte Carlo simulations. In this section we compare our results with results of MC calculations obtained by 23 and 24 
for a single polymer chain trapped inside a slit of two repulsive walls, what corresponds to a canonical ensemble. The 
canonical free energy can be obtained via a Legendre transform from the grand canonical one in the thermodynamical 
limit (N, V — > oo) in the form: 

F{Ni) = n\f i {Ni)]+ t i{N l )Ni, (7.1) 

with Q from (3.9). 

Thus the reduced canonical force for a one polymer chain Ni = 1 can be written as its dimcnsionless counterpart: 

= — tt(Luj). (7.2) 

k b T udL y ' K ' 

It should be mentioned, that both Monte Carlo algorithms (see 23 and 24 ) differ very much from each other in the 
range of analyzed slit widths and chain lengths of the simulated polymers. In 23 an off-lattice bead and spring model 
for the self-avoiding polymer chain in d = 3 dimensions trapped between two parallel repulsive walls at distance D 
has been studied by Monte Carlo methods, using chain lengths up to N < 512 (number of monomers in the chain) 
and distances D from 4 to 32 (in units of the maximum spring extension) . It was stated that the total force K exerted 
on the walls is repulsive and diverges for the case of narrow slit as 

S? -<£>-""■ ™ 

where R g is the radius of gyration of the polymer chain in unrestricted geometry. 

In Ref. 24 a lattice Monte Carlo algorithm on a regular cubic lattice in three dimensions, with D lattice units in 
z-direction and impenetrable boundaries was applied. The other directions obeyed periodic boundary conditions. 
The proposed MC simulations 24 are based on the analytical result obtained by 17 for the scaling behavior of partition 
function for a chain confined to a slit geometry of width D: 

Z N (D) oc ^ oc + aD-^y N N^- 1 D^ 22 ^ 1 , (7.4) 

for N, D — > oo, but D <C N v , where (j,^ is the critical fugacity per monomer and 7^ is the universal exponent (see 
(2.1)) dependent on the space dimension d and the parameter a is a universal amplitude. The critical fugacity means 
the averaged inverse number of possible steps at each site. In 24 universal amplitudes and exponents for the partition 
function of a chain trapped in the slit with respect to that of a free chain have been obtained through analyzing the 
statistics for different D and number of chain monomers up to N < 80 000. Also both cases of ideal chains (modeled as 
a simple random walk (RW)) and chains with excluded volume interactions (modeled via self-avoiding walks (SAW)) 
have been studied. In the case of an ordinary random walk on a regular cubic lattice in three dimensions one has 
obviously fi^ — g and 7^=3 = 1. In the case of a SAW on such a lattice it is clear that at least fi^ > -|. From Eq. 
(7.4) one may obtain the force exerted onto the walls in units of ks T as: 

K =k B T-^ In (Z N (£>)). (7.5) 

In the limit (D < N v , N,D — > 00) K becomes: 

K = k B T(^D-^) . (7.6) 

One should note that all functions here are in terms of dimensionless length scales, the number of lattice sites (D and 
N). In order to compare with our results it must be translated into terms of L and R g . Apparently L = uD, with u 
denoting the lattice spacing, and the reduced, dimensionless force reads: 



KD a /_6_V" 3 (L_ 



k = ^ = -^-[- 2 ) ' (^) , (7-7) 
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where we take into account the general relation (e.g. 25 ): 

9 ob 2 N 2v * , N 

R g = xi— ^— (7.8) 

in d = 3 dimensions. Parameter b denotes the (effective) segment length of the polymer model under consideration. 
In the case of RW and SAW on a cubic lattice one has simply b — u because the segment length in these models is 
naturally provided by the lattice parameter u. It? 4 the universal amplitude a for the case of ideal chains was found 



Comparison with MC simulations - 2 repulsive walls 
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FIG. 5: Comparison of our theoretical results with Monte Carlo simulations for a trapped chain between two repulsive walls 

The plots Ideal chain (exact) and EVI (wide slit) represent the results of our calculations. MC - RW and MC - SAW are due 
to the estimated asymptotic behaviour in the narrow slit limit by 24 for random walks and self avoiding walks. MC - EVI are 

the results obtained by 23 . 

as a ~ 0.2741, which is very close to the exact value, computed analytically in 17 , of a — ^g. Taking into account that 
for ideal chain v = 0.5, Xd — 1 and /loo = | the force becomes: 

k ld = 2tt 2 (J^j . (7.9) 

In Figure 5 this asymptotic behaviour for narrow slits is clearly recovered by our results for ideal chains, where the 
narrow slit limit is valid. By contrast, for a SAW in Ref. 24 was suggested the value a « 0.448. Taking into account 
the values for v s=a 0.588, %3 ~ 0.958 25 and w 0.2135 the reduced force can be written as 

L ^ - 1 - 7 



k saw « 16.95 ^— J . (7.10) 

The result Eq.(7.10) is presented in Figure 5 in its regime of validity and compared to our theoretical results for a 
trapped chain with EVI, which are valid for the wide slit regime. As it easily can be see from Figure 5, the result 
(7.10) very well fit to our predictions in wide slit limit. Also, in Figure 5 the results obtained by the authors of Ref. 23 
are plotted and one notes a qualitative agreement to our predictions. One of the possible reason for the remaining 
deviations with results of Ref. 23 is that the chain in the MC simulation is too short in order to compare with results 
of field-theoretical calculations. It should be noted, that at the moment no simulations concerning two inert walls or 
one inert/ one repulsive wall exist. 



C. Experiment 

In Ref. 8 an experimental study of the depletion effect between a spherical colloidal particle immersed in a dilute 
solution of nonionic linear polymer chains and a wall of the container through total internal reflection microscopy was 
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analyzed. Using the Derjaguin 11 approximation we could compare our theoretical results with experimental data in 
the case when the radius of the spherical colloid particle R is much larger than radius of gyration R g and the closest 
distance a between particle and the surface. The deviation of the experimental setup to the presented theoretical 
approach connected with the fact that the second wall is not plane but curved. Summing up the depletion potential 
per volume unit for the case of two plane surfaces in the margins of the curved volume allows to estimate the depletion 
effects in the case of a sphere and a wall. In the experiment by 8 the radius of gyration was measured as R g = 0.101 \im 



and the colloidal particle was reported to have a radius R 
approximation yields: 



1.5 fim. Straightforward application of the Derjaguin 1 



4>depl (a) 
n p k b T 



= 2irRi 



dy (R 



Rxv) e(y) 



(7.11) 



with a the minimal distance between the sphere and the wall. Since in the range of y the last two terms in the 
parenthesis are much smaller in comparison to the first one we can assume that: 



4>depl{a) 



2nRR: 



oo 

J dye(y) . 



(7.12) 



n p k b T 

The experimental data in comparison to our theoretical prediction are plotted in Figure 6. It should be mentioned, 



E 
i 





-0.05 
-0.1 

-0.15 
-0.2 

-0.25 
-0.3 



Comparison with * dep | from [RBL98] 



EVI with R =0.101nm 
experimental data [RBL98] 
Ideal chains with R„=0.101nm 
EVI by [SHKD01] 



0.05 



0.1 0.15 
a [urn] 



0.2 



0.25 



FIG. 6: Comparison of approximated theoretical results with experimental observation, due to 8 



that our results obtained in the framework of the massive field theory are situated slightly closer to the experimental 
data than previous theoretical results obtained in the framework of the dimcnsionally regularized continuum version 
of the field theory with minimal subtraction of poles in e = 4 — d ls . Unfortunately, this shift is not enough in order 
to obtain quantitative agreement with experimental data. But, the obtained theoretical curves in Figure 6 are in 
qualitative agreement with experimental data. The quantitative discrepancy can be removed if we use the radius of 
gyration as adjusting parameter by analogy as it was done in 18 . From another side this indicate about importance 
of further theoretical investigations of depletion interaction potential and depletion force in the crossover region from 
wide to narrow slit. 



VIII. CONCLUSIONS 



Using the massive field theory approach directly at fixed dimensions d = 3 we calculated the depiction interaction 
potential and depletion force between two repulsive, two inert and one repulsive and one inert walls confining a dilute 
solution of long flexible polymer chains. The obtained calculations for all cases of polymer-surface interactions were 
performed for the ideal chain and real polymer chain with excluded volume interactions in the wide slit regime. 
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Besides, we used some assumptions which allowed us to estimate the depletion interaction potential in the region of 
narrow slit. Our results are obtained up to the next e~ 2fJ,L order in comparison with results of e-cxpansion 18 . Our 
investigations include modification of renormalization scheme for the case of two inert walls (or mixed walls) situated 
on big, but finite distance L with L > R g such that the polymer chain is still not deformed too much from its original 
size in the bulk. The obtained results indicate that the reduction of the depletion effect due to excluded volume 
interactions is less stronger within the massive field theory approach as compared to the dimensionally regularized 
continuum version of the field theory with minimal subtraction of poles in e = 4 — d 18 in one- loop order. We found very 
good agreement with Monte Carlo simulation data 24 and 23 for the case of two repulsive walls. Taking into account 
Derjaguin approximation we obtained good qualitative agreement with experimental data 8 for the depletion potential 
between a spherical colloidal particle of big radius and repulsive wall. From comparison of obtained theoretical results 
and experimental data we can see that the results obtained in the framework of the massive field theory are situated 
slightly closer to experimental data. But, this shift is not enough in order to obtain good quantitative agreement 
with experiment. Interesting fact is that even the taking into account the excluded volume interaction between the 
monomers of the polymer chain do not resolve completely this problem. One of the possible ways to find a good 
agreement could be connected with further theoretical investigation of crossover region from wide to narrow slit. 
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Appendix A: The surface contributions 

To calculate the function T defined in Eq.(3.14), we need the free propagator for a semi- infinite system confined by 
a surface at z = 0. This free full propagator has a form 21 : 

G* i:jH s(p,p',z,z';^o,co) = (2n) d - 1 5 ll 8(p + p')G H s{p, z, z! , Mo, c ) , (Al.l) 

with 

G HS (p, z, z', Mo , co) : - ^- (e-^ z - z '\ + ?Ll* e -«°{*+>')\ ; (AL2) 



where k = \J p 2 + 

In the zero-loop order we have: 



T l = \^. (A1.3) 

In one-loop order the calculation for Dirichlet boundary conditions on the surface (or ^ — > oo) yields after renormal- 
ization in fixed dimensions d — 3: 

r D = ±(l-^vlnl). (A1.4) 
yu 3 \ n + 8 8 / v ; 

And for Neumann boundary conditions (c = 0) after renormalization we obtain: 

--4(ln2-i)^± 2 
M V 2) n + 



T In2-- (A1.5) 



where we introduced rescaled renormalized coupling constant v in the form: v = ^"g 8 ^ ff^J/l ^- The correspondent 
fixed point in one- loop order approximation is v* =1. 
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Appendix B: The case of finite slit separation for two inert walls 

Taking into account the new <5cf" s (see Eq.(4.9) we can calculate Sf in accordance with Eq.(6.12) for the case of 
big, but finite slit separation L. We obtain: 



2 A< 5 ~ 5 ) c- 2 ^ 



<f/« XC , 2 ^ ( -j + ) . (A2.1) 



After substitution of A (s s) from (4.10) the result for Sf is: 

Sf w - XC^ b| | + C E + 7 In 2 - 6 + In M L - e 4 ^ L Ei(-4 M L)^ } ■ ( A2 - 2 ) 

If we carry out the inverse Laplace transform, the result for <d(y) in the wide slit limit is: 



90/)~- Cg + 7 2 ln2 ~ 7 (/f e " V - 2 y orfc(V2^ - ^crfc(V2 y ) 



- ! + f (^ Ei (-4V7)j . (A2.3) 

In contrast to (6.11) this expression is indeed negative. Thus, if we perform calculations for the depletion interaction 
potential and the depletion force including big, but finite slit separation L we obtain, that force for the case of two 
inert walls change character and becomes attractive. In Figure 3 the depletion interaction potential and the depletion 
force obtained in the framework of this alternative renormalization scheme with Sc s m are plotted in comparison to 
the results obtained via the original renormalization using Sc sem i. Here we introduced for convenience the following 
notations: Sc sem i = Sci and Scf ~ s = Sc s u t with i = 1,2. 
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